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Abstract 

We study the dynamics of renormalization of a specific interval exchange trans- 
formation which features exact scaling (the cubic Arnoux-Yoccoz model). Using a 
symbolic space that describes both dynamics and scaling, we characterize the peri- 
odic points of the scaling map in terms of generalized decimal expansions, where the 
base is the reciprocal of a Pisot number and the digits are algebraic integers. The set 
of periodic points has a rich arithmetic and geometric structure: we establish rigor- 
ously some basic facts, and use extensive numerical experimentation to formulate a 
conjecture. 
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1 Introduction 



There are symbolic representations of dynamical systems which lead to interesting arith- 
metical phenomena. In such arithmetic codings, the symbols represent generalized 'digits' 
in a number system which naturally describes the dynamics. For a survey of this rich area 
of research from a dynamical perspective, see [32]; an arithmetical viewpoint is developed 
in [1-3]. 



In arithmetic codings — even in the simplest settings — one encounters at once inter- 
esting problems. A well-known example is the doubling map of the unit interval (circle) 
7 : x i— > 2x (mod 1), with the symbolic dynamics defined by the partition of the interval 
into two equal halves. The symbols representing an orbit are just the binary digits of its 
initial point, so that the eventually periodic symbolic sequences correspond to the rational 
numbers, and the periodic sequences to the rationals with odd denominator. For each n, 
the periodic set 

Fix ^ n ) = ^T Zn[0 ' 1) n ~ 1 (1) 
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is a lattice. The exponentially large denominator determines the maximal complexity 
(height) an n-cycle can have. 

The question of minimal complexity is altogether more difficult: when 2 n — 1 has 
a large number of divisors, cancellation becomes likely, leading to n-cycles with much 
smaller denominator. This phenomenon becomes evident when one looks at the twin 
asymptotic problem, namely that of determining the period of rational points with a given 
odd denominator. In this setting, as the denominator increases, the low-complexity cycles 
are those that are found first. It is possible to show that this problem is closely related 
to the so-called Artin's conjecture on primitive roots [29]; at present, the asymptotic 
behaviour of the period of a rational point can only be established assuming the validity 
of the generalized Riemann hypothesis [25]. 

In the study of periodic orbits, a natural generalization of the above construction 
occurs in hyperbolic toral automorphisms. In place of Q, we now have the field Q(A), 
where A is an eigenvalue of the automorphism, while the lattice ([1]), which is a Z-modulqj, 
is replaced by a Z [A] -module I n . The dynamics is given by modular multiplication by A, 
and the study of periodicity involves asymptotic problems that are quite similar to those of 
the rational case [6,13,27]. The symbolic representation of hyperbolic toral automorphism 
was considered in [17], and is related to the codings corresponding to expansions to non- 
integral bases, the so-called beta- expansions [2,3,32]. We note that a module structure for 
periodic points is also found for certain dynamical systems over finite fields. The set of n- 
cycles of a linearized polynomial is a vector space, that is, a module over a field [5, lemma 
4.2] (see also [12]). 

Of a different nature are the codings related to renormalization, which first appeared 
in dynamics in the study of persistence of planar rotational motion under perturbation 
[23,31]. In this case the scaling process is a second-order recursion over a suitable space 
of maps, and the symbols associated to it are the coefficients of the continued fraction 
expansion of the rotational angle. As a result, the renormalization dynamics is eventually 
periodic when the rotational angle is a quadratic irrational. (These dynamical systems 
should not be confused with the so-called rotational expansions [32].) 

In this paper we consider a specific one-dimensional dynamical system — an interval 
exchange transformation proposed by Arnoux and Yoccoz — which features exact scaling, 
meaning that the induced map (Poincare section) on one of the intervals is a scaled version 
of the original map. This map is defined over a cubic number field — see equation ([3]). An 
interesting feature of this model is the existence of a symbol space which can be used to 
describe both the dynamics of the map and the dynamics of scaling. The dynamics of the 
map has zero topological entropy, and its symbolic representation could be described as 
a generalized odometer [22]; by contrast, the dynamics of scaling is a shift with positive 
entropy. 

1 given a ring R, an 7?-module is an an additive group equipped with multiplication by elements of R 
— see [15] 

2 A polynomial over a finite field is linearized if all exponents are powers of the characteristic of the 
field. 
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This type of coding originates from the so-called recursive tiling property of a dynamical 
system: a scale-invariant set of tiles exists, which cover the entire space under iteration, 
in a hierarchical fashion. The symbol space consists of pairs of integers, which label a tile 
together with its transit time measured from a reference position. These tilings have been 
used extensively in the study of piecewise isometric systems in two dimensions [18,20,21]. 
It turns out that these symbol sequences posses a rich arithmetical structure, which, to 
our knowledge, is largely unexplored. In this respect, the Arnoux-Yoccoz model is an ideal 
object of study; it illustrates how basic constructs developed elsewhere generalize, while 
remaining sufficiently simple for the numerical exploration of its asymptotics. 

This paper is devoted to the study of the scaling dynamics — a subshift of finite type — 
and in particular of its periodic orbits. This mapping leads to an arithmetic coding where 
the digits are algebraic integers rather than integers (a phenomenon already studied in 
unrelated contexts [14,17]) and where the eventually periodic sequences coincide with the 
elements of the underlying field. 

The system under study is defined in section [2j where some key results in [22] are 
reviewed. In section [3] we show that the periodic points belong to the cartesian product 
of a module over a finite ring, which depends on the period, and a finite set — the core 
region, which is period-independent. These two components represent the 'fractional and 
integer parts' of the periodic points, respectively. The finiteness of the core region derives 
from the Pisot property of the scaling factor. In section HI we order the periodic points by 
increasing period, and consider asymptotic properties of their integer and fractional parts. 
We prove that the set of periodic points which share the fractional part with some other 
periodic point is infinite. Numerical evidence, however, strongly suggests that this set has 
zero density, which is expressed as conjecture [TJ By embedding the fractional parts of the 
cycles into the unit cube, we study the relation between fractional and integer parts from 
a geometrical viewpoint. We collect some evidence that the phase space is organized in a 
highly nontrivial manner, which suggests directions for future investigations. 

In section [5j we order the periodic points by increasing 'denominator', and relate their 
period to the ideal factorization of the denominator in the ring of integers of the underlying 
number field (proposition H]). Our experiments here investigate the asymptotic distribution 
of periods for large denominator, obtaining a picture consistent with that of increasing 
period. 

This research was supported by EPSRC grant No GR/S62802/01. 



2 Preliminaries 

Arnoux and Yoccoz [4], in their study of pseudo-Anosov diffeomorphisms, introduced a 
family of interval-exchange transformations defined over algebraic number fields of arbi- 
trarily large degree, with the property of being renormalizable, namely of generating only 
finitely many induced maps, up to scaling. The simplest nontrivial member of this family 



3 



corresponds to a cubic field; the object of our study is the scale-invariant (exactly renor- 
malizable) mapping p induced by it on a sub-interval. This mapping p : [0, 1) — * [0, 1) is 
defined as follows [4, 22] 

p(x) = x + Tj, x £ Qj VLj = [5j, dj+i) j = 0, . . . ,N — 1 (2) 

where N = 7, the discontinuity points Sjy = 1 and translations <L-, Tj are listed in table [TJ 
and A is the real root of the irreducible cubic polynomial 

f(x) = x 3 + x 2 + x - 1. (3) 

The map p corresponds to the permutation (0,1,2,3,4,5,6) i— ► (6,2,4,3,5,1,0) of the 



j 




T 3 


v i 


P(j, t) 








A + A 2 


4 


(0,6,3,6) 


1 


1 - A - A 2 


-1 + 3A 


13 


(0,6,3,6,1,6,2,5,6, 1,6,3,6) 


2 


1 - 2A + A 2 


A -A 2 


12 


(0,6,3,6,1,6,2,5,6,2,4,6) 


3 


(3 - 4A - A 2 )/2 


-1 + 2A + A 2 


8 


(0,6,3,6,1,6,3,6) 


4 


(1 - A 2 )/2 


1 - A - A 2 


8 


(0,6,4,5,6,2,4,6) 


5 


(-l + 2A + 3A 2 )/2 


A- A 2 


12 


(0,6,4,5,6,2,5,6,1,6,3,6) 


6 


A 


-A 


4 


(0,6,4,6) 



Table 1: The data defining the cubic scale-invariant Arnoux-Yoccoz map p. 

N disjoint intervals Oj, which constitute a partition of [0,1). The resulting dynamical 
system has zero topological entropy, being a piecewise isometry [9]; in addition, it is 
uniquely ergodic, and hence, in particular, it has no periodic orbits. 

The map p is defined over the cubic number field 

Q(A) = {r + r 1 X + r 2 X 2 \r l E Q}. (4) 

If we restrict the coefficients in (JH) to integer values, we obtain the ring Z[A] of all 
algebraic integers in Q(A); its group of units (invertible elements) has rank one, and is 
generated by A [11, p. 519]. The unit u = A 3 is important: in fact, scaling by uj conjugates 
the map p to the induced map on the leftmost interval £1$ = [0, 1 — A — A 2 ) = [0, to), which 
can be verified by direct calculation. 

To the above data we associate the sequence space V constituted by the sequences 
a = ((j 1 ,t 1 ),(j 2 ,t 2 ),...) 0<j k <N, 0<t k <v k 

where the jk are subject to the constraints 

jk=P(jk+i,tk+i), fc >! (5) 

and p is the path function defined in table [TJ Each sequence a £ V satisfying © defines 
a real number x via the equations 

oo ti—1 

i=l t=0 
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where the 'digit' di = di(ji,ti) belongs to a finite subset of the Z-module M generated 
by the translations tj-. It is easy to see that Ai = Z[A]; furthermore, there are ^ uj = 61 
valid symbols, corresponding to 25 distinct values of the digits. 

By construction, the point x(o~) belongs to the unit interval. The set of points sharing 
a finite digit sequence (di,...,d n ) is a half-open interval, called a tile of level n: it is 
the intersection of translated and scaled versions of the original intervals Qj. For each 
integer n, the tiles of level n constitute a partition of the unit interval: this is the recursive 
tiling property of the mapping, which implies the surjectivity of the function a \— ► x{a). 
The representation ([6]) may be regarded as a variant of the expansion of real numbers in 
non-integral bases — the so called /3-expansions. Our basis u is the reciprocal of a Pisot 
numbeio, a property which is of importance here (it underpins theorems HJ El below) as 
it is in the theory of /5-expansions [1,32]. We emphasize again that the digits of our 
w-expansion are algebraic integers. 

The map a i— > x(a) is not injective, but the lack of injectivity is easily controlled. 
Whenever it fails, there are precisely two codes corresponding to the same point, and 
these codes are eventually periodic and entirely characterized by their periodic part. This 
degeneracy can be removed by stipulating that the nested set of tiles specified by the code 
a must also contain its limit point x. If only such codes are allowed, the correspondence 
between points and codes is bi-unique. This amounts to excluding the following six periodic 
tails 

(1,9)°° (2,6)°° (3,2)°° (4,6)°° (5,3)°° (6,1)°°. (7) 

(The situation is analogous to the decimal representation of rational numbers: the tail 
(9)°° is excluded to achieve a unique representation.) 

With the above stipulation on valid codes, the map 

7 : [0,1) -» [0,1) x^{x-di)uj- 1 (8) 

is conjugate to the left shift on V. This is a subshift of finite type, which we call the 
scaling dynamics of the Arnoux-Yoccoz map. The transitivity of the incidence matrix 
corresponding to the admissibility conditions ([5|) can be verified directly. It then follows 
[16, p. 51], that 7 is mixing, and that the periodic orbits are dense. 

The following result was proved in [22]. 

Theorem 1 A point x of the unit interval belongs to Q(A) if and only if the code o~(x) is 
eventually periodic. 

This result calls to mind the familiar property of expansions in integral bases, where the 
set of points characterized by the eventual periodicity of the digits is the field Q. To place 
this theorem into context, we note the general results of Bertrand [7] and Schmidt [30] on 
eventually periodic greedy expansions in a Pisot base. 

3 A Pisot number is an algebraic number greater than one, whose algebraic conjugates other than itself 
lie inside the complex unit circle. 
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To study dynamics over Q(A) we consider the decomposition 

Q(A) = E + Z[A] 

where 

s = {eo+6A + e 2 A 2 ieiGQn[o,i)}. 

This prescription gives a unique representation of x G Q(A) as x = £ + /?, with ^GS and 
/? G Z[A]. Accordingly, we define the projections to the two components 

11/ : Q(A) -» E x^£, n, : Q(A) -> B x ^ (3. (9) 

These components should be viewed, respectively, as the 'fractional and integer parts' of 
the field elements. Note that S is a Z[A]-module, isomorphic to Q(A)/Z[A]; this fact will 
be of importance later. 

In the rest of this paper we shall be concerned with the set of periodic points of 7, a 
7-invariant set which belongs to Q(A) n [0, 1). 



3 Periodic points of scaling dynamics 

Consider a point x with strictly periodic code 

<r(x) = ((ji,h),...,(j n ,t n )r. 

From theorem [T] we have that x 6 Q(A) D [0, 1), and a straightforward calculation shows 
that 

1 n 

i=l 

where the digits dj were defined in §6§. The following result — proved in [22] — is crucial 
to our investigation 

Theorem 2 There exists a finite subset B o/Z[A] such that any ^-periodic point x can be 
represented as 

x = i + (3 CeS,/?eB. (11) 

Conversely, for every £ e H, there exists ft E B such that x = £ + (3 is periodic under 7. 

The set B G Z[A] is called the core region. The theorem says that the points which are 
periodic under the scaling map 7 can have only finitely many distinct integer parts. Exact 
computations give the estimates 

31 < #B < 348. (12) 

The upper bound was established using a technique developed in [22] . (In that reference, an 
upper bound of 225 was obtained for the projection of the core region to a two-dimensional 
lattice.) The lower bound will be established below — see table [3j 



6 



Theorem [2] outlines the basic structure of the periodic points of 7, yet important ques- 
tions concerning their integer and fractional parts are not addressed. Are there infinitely 
many fractional parts £ G S shared by more than one periodic point? Is it true that for 
all (3 in the core region there are infinitely many periodic points with integer part equal 
to some (31 Furthermore, the above questions can be refined to formulate statements on 
the density of the relevant infinite sets. 

To investigate these issues, for (3 G B we consider the set of all £ G E for which 
£ + (3 is periodic. Using the projections defined in (jSJ), we see that 

H(/3) = n / on^ 1 (/3) (13) 

and from theorem [2] we know that 

U = S - ( 14 ) 

/3eB 

Next, for £ G E we define the multiplicity to be the number of 7-periodic points with 
fractional part equal to £. From theorem [2] we know that 1 < < and we are 
interested in the set <9E of those £ that correspond to more than one periodic point. This 
set is given by 

dZ = {£ e H : «(0 > 1} = |J (H(/3) n B(^)) . (15) 
For reason that will become clear later, we call dS the boundary set. 



We begin our study of periodic points by seeking an analogue of equation ((TJ. Let M n 
be the smallest positive integer which is divisible by 1 — uj u in Z[A], that is, 

{l-u n )a n = M n (16) 

for M n minimal and some a n G Z[A]. From equation (fTOl) one sees that M n x G Z[A], that 
is, 

Equivalently, M n is the least common multiple of the denominators of all n-periodic points. 
In practise, M n is computed efficiently with the extended Euclid's algorithm. 



By construction, M n is a divisor of the norm N(l — u; n ) of 1 — iv n (the norm of an 
algebraic number is the product of its algebraic conjugates). While the value of N(l — 
uj u ), is well-behaved — see section U equations (I20|2ip — the value of M n features large 
fluctuations. In this respect, we note the factorization 

l-u n = -HC 3d (X)l[C d (X) (17) 

d\n d W 
S/d 

where Cd(x) is the d-th cyclotomic polynomial [26, p 37]. This formula is established 
starting from the polynomial identity x n — l = Y\ d \ n C d {x), recalling that uj = A 3 , and then 
relating the divisors of 3n to those of n. The formula provides only a partial factorization 
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of 1 — ui n , because the irreducibility of the cyclotomic polynomial Cd(x) in Z[x] does not 
imply that Crf(A) are primes in Z[A] (which is a principal ideal domain [11, p. 519]). 

For example, the identities 

u - 1 = Ci(A)C 3 (A) = (A - l)^ 1 (1 - A) 3 A~ 5 = 2 

together with the fact that A is a unit, show that the denominator of a periodic point 
always has the irreducible factor 1 — A, which is a prime divisor of 2. From the above 
observation, it follows that M n is even. Likewise, Cq(X) is a divisor of 7, and so if n is 
even, M n is divisible by 14, etc. 

The periodic orbits of the scaling dynamics can be ordered in two ways, namely by 
increasing period n, or by increasing denominator m (the choice of sub-orderings for given 
n or m are not important here). As mentioned above, these orderings lead to different 
perspectives (see also [10, section 6]), and they will be considered in sections 0] and 
respectively. In analogy with the rational case, whenever 1 — of 1 has a large number of 
divisors, cancellation becomes likely in equation (|10p . leading to n-cycles with denominator 
much smaller than M n . As the denominator is increased, these are the cycles that are 
found first. 

Each ordering gives asymptotic information about the sets H(/3) and dE defined above. 
This will be our main object of study. 



4 Ordering by period 



Let the period n be fixed. With reference to equations ([9]) and (|TTj) . we define 

Zn = % (—L^Z[\]\ T' n = U f (Fix( 7 «)) B n = IT, (Fix( 7 n )) . (18) 

From equation (|10p it follows that the fractional part of the n-cycles is contained in 2 n , 
and hence X' n Cl n . By construction, the fractional parts in the (possibly empty) residual 
set T n \ X' n belong to cycles whose period is a multiple of n. 

To determine the cardinality of Z n , we must refine the decomposition (llip . Let m be 
a positive integer, and let H m be the set of ^-points with denominator m namely 



(19) 



which is again a Z[A]-module, and is isomorphic (via the map : x i— > mx) to the module 
Z[A]/mZ[A], which has m 3 elements. 

Consider now the algebraic integer a n defined in equation (|16p . The map ip sends I n 
to the ideal generated by a n in Z[A]/M n Z[A] (now regarded as a ring), which yields the 
equation 

#E Mn = #l n x\N(a n )\. 
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Using the multiplicativity of the norm, and the fact that #H M = N(M) = M 3 , we obtain 
the formula 

#X n = \N(l-u; n )\ (20) 

which, together with equation (I16p . shows that #2T n is a multiple of M n . Explicit compu- 
tation gives 

#X n = uj- n - uj n - 2 - v 7 ^™) cos(n#) (21) 

where 

cos(#) = -\/^ (5 — w). 
To compute $ Fix(7 n ) we consider the incidence matrix 



/ 1 








1 








2 \ 


1 


2 


1 


2 





1 


6 


1 


1 


2 


1 


1 


1 


5 


1 


1 





2 








4 


1 





1 





2 


1 


3 


1 


1 


1 


1 


1 


2 


5 


V i 











1 





2 / 



where Ajj is the number of times that the first-return orbit of the scaled interval wOj 
visits the interval Oj, obtained from the path function data of table [TJ The characteristic 
polynomial of A factors into irreducibles as (x — l)(x 3 — 5x 2 + 7x — l)(x 3 — 7x 2 + 5x — 1), 
whose real roots are 1, u> = A 3 , and uj -1 , respectively. Accounting for the six forbidden 
period-1 codes listed in 0, we find 

# Fix( 7 n ) = TrA n - 6 = uj~ n + to n + 2(v^F™ + y/uF) cos(n6») - 5. (22) 

Expressions (|2ip and (I22p have the same leading term uj~ n (recall that \u\ < 1), and hence 
# Fix(7 n ) ~ #X n . Furthermore 

# Fix( 7 n ) - #J n = 4\/^ r " cos(n6») + 2w n - 5 ~ 4\/^F™ cos(n6>). (23) 

This formula shows that for some values of n, ^Fix(7™) is greater than which 
implies that there exist distinct periodic points with the same fractional part £. In fact, 
the asymptotic expression (|23p shows that this must happen infinitely often, and that the 
population of points with multiplicity greater than 1 is maximal when the period is a 
denominator of the continued fraction expansions of 9/2n. Recalling the definition (|15p of 
the boundary set, we have the following result. 

Proposition 3 The boundary set <9H is infinite. 

Using induction, it is possible to prove that the periodic codes 

((1, A)\ (3, 4), (1, 2)) 00 ((2, 9) fc , (4, 5), (2, 10))°° k = 0, 1, . . . 

(for k = the symbol is deleted) correspond to points of period + 2 on the boundary set. 
Now, the fixed points (1, 4)°° and (2, 9)°° also belong to dE, as easily verified. This shows 
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that the boundary set contains points of any period, which is stronger than proposition [3j 
We shall not produce this proof here. 



When #Fix(7™) < #T n , the ^-values do not exhaust the whole of T n , and formula (|23|) 
again shows that this must happen for infinitely often. Since ^I' n ^B > #Fix(j n ), we 
obtain 

This bound is far from optimal. To see this, we define 

Constructing explicitly all n-periodic points for n < 14, gives the data displayed in table 

m 



n 




#X" n /# Fix( 7 n ) 




1 


2 


0.2857 


7 


2 


14 


0.3111 


17 


3 


38 


0.1467 


21 


4 


86 


0.05595 


22 


5 


182 


0.01939 


24 


6 


374 


0.006440 


24 


7 


758 


0.002100 


27 


8 


1534 


0.0006831 


30 


9 


3170 


0.0002268 


30 


10 






30 


11 






30 


12 






30 


13 






30 


14 






31 



Table 2: Multiplicity and core region data from orbits of period n. 



The multiplicity data are limited to the range n < 9, because the comparison of all 
fractional parts requires that these quantities be stored. The most significant finding, to be 
used in conjunction with the asymptotic formulae (I2ip and (I22p . is that distinct periodic 
points typically have distinct fractional parts; namely, the set of cycles with multiplicity 
1 has full density. Furthermore, the ratio ^X" n /# Fix(7 n ) appears to be tending rather 
rapidly to zero, with a scaling ratio close to 1/3. In the computations to be described 
in section [5] we observe the same phenomenon when the cycles are ordered by increasing 
denominator, leading to the following conjecture. 

Conjecture 1 The following holds 

Hm #{g£X' n : «(Q = 1} = nm #{g £ E m : «(Q = 1} = ^ 

rwoo m-»oo m? 

The validity of this conjecture would imply that, as n — > oo 
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to be compared with the bound ([24"]) . 

Regarding the core region data of table EJ the lower bound of 31 is substantially smaller 
than the upper bound of 348 given in (|12p . With reference to equation (jlip . we display 
in table [3] some elements = mo + mi A + 7712A 2 of the core region, represented as integer 
triples (mo, mi, m2). Each value of (5 has a probability defined as the limiting density 
of the ^-values that correspond to it. Approximate values of these probabilities may be 
computed in two different ways, corresponding to ordering by period (the second and third 
columns in the table) and by denominator (the fourth column). In the former case, the 
data was computed from the 14-cycles (130399019341 data points), which project to the 
31 points of the core region listed in the first column, respectively. The data for period 
13 are also shown, to give an idea of convergence. The probability fj,' computed using 
all cycles with 'denominator' m < 200 is displayed in the third column (373112717 data 
points, which project to a subset of 28 core region points) — see the next section. While 
it is conceivable that the core region could contain more points than those displayed in 
table [3l the stability of the above figures suggests that they are a reliable approximation 
to the densities. 

So far we have expressed probabilistic information in terms of densities; we close this 
section by discussing related questions of measure. To this end, we consider the embedding 
of the fractional parts into the unit cube 

ip:Z^{0, if r + r 1 X + r 2 X 2 ■-> (r ,ri,r 2 ), < n < 1. (25) 

This construct depends on the particular basis chosen for Q(A). However, a change of 
basis corresponds to a unimodular transformation, which, if we identify the unit cube 
with the 3-torus, is continuous and volume preserving. So any property of the (closure of) 
embedded objects which is topological or concerns three-dimensional Lebesque measure 
will be independent of the basis. 

It is natural to consider sequences of points that are uniformly distributed in the cube; 
for instance, the sequence of fractional parts of periodic points with increasing period, or 
with increasing denominator. Consider now the embedding y?(9H) of the boundary set 
<9E, defined in equation (fT5|) . Recalling that #Fix(7 ra ) ~ uj^ 1 = A~ 3 (see equation (|2"2"|) ). 
the scaling law observed in table [2] is consistent with the closure of <p(dS) having zero 
three-dimensional Lebesque measure. Now, the points of the boundary set are precisely 
the points with multiplicity greater than one; so any uniformly distributed sequence will 
have the property that the density of the elements with multiplicity greater than 1 is zero. 
Thus the vanishing of the measure of the boundary set in the cube implies conjectured! 

We have examined the arrangements of the boundary set in the cube, which gives some 
support to this conjecture, although the available data are not quite conclusive. In fact, 
there is an even stronger property that should be considered, namely that the core region 
localizes in the unit cube; this means that the measure of the closure of the sets <p(5(/3)) 
add up to unity, and represent the probabilities /u(/3) estimated in table [3l (Loocalization 
techniques have proved very useful in the study of certain lattice maps [8, 19].) 

These problems deserve further investigation. 
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p 




A*13 




M14 




A* 




( o, 


o, 


o) 


0.21275 




0.21275 




0.21566 




( o, 


-1, 


-1) 


0.18921 




0.18921 




0.18811 




( o, 


o, 


-1) 


0.13777 




0.13777 




0.13688 




( 1, 


-2, 


-1) 


0.095625 




0.095625 




0.096146 




(-1, 


o, 


o) 


0.077156 




0.077155 




0.075485 




(-1, 


o, 


1) 


0.062212 




0.062213 




0.062157 




( o, 


-2, 


0) 


0.57953 




0.057952 




0.057331 




( o, 


-1, 


1) 


0.029309 




0.029308 




0.029370 




( 1. 


-3, 


0) 


0.026513 




0.026514 




0.026326 




( 1. 


-2, 


0) 


0.026416 




0.026416 




0.027329 




( o, 


-2, 


-1) 


0.022015 




0.022014 




0.020936 




( o, 


-2, 


1) 


0.013304 




0.013304 




0.013345 




( o, 


-1, 


o) 


0.011008 




0.011009 




0.010908 




( o, 


o, 


1) 


0.0088052 




0.0088053 




0.0093162 




(-1, 


1, 


o) 


0.0079820 




0.0079830 




0.0083504 




(-1, 


1, 


1) 


0.0068145 




0.0068150 




0.0069532 




( o, 


1, 


o) 


0.0056184 




0.0056188 




0.0058062 




( 1. 


-1, 


-1) 


0.0052609 




0.0052609 




0.0056897 




( 1. 


-3, 


-1) 


0.0020413 




0.0020418 




0.0019191 




(-1, 


o, 


-1) 


0.0011704 




0.0011707 




0.0010408 




(-1, 


-1, 


1) 


0.00075012 




0.00075025 




0.00065573 




( 1. 


-4, 


o) 


0.00020058 




0.00020069 




0.00018132 




( o, 


-3, 


o) 


7.5098 • 10" 


-5 


7.5092 • 10" 


-5 


5.2630 • 10- 


-5 


( 1. 


-2, 


1) 


2.1778 • 1(T 


•5 


2.1757 • 10- 


-5 


3.2518 • 10- 


•5 


( 1. 


-3, 


1) 


3.7165 • 10" 


-6 


3.7439 • 10" 


-6 


7.8207 • 10- 


-6 


( o, 


-1, 


-2) 


3.6249 • 10" 


-6 


3.6509 • 10- 


-6 


1.1444 • lO - 


-6 


( 1. 


o, 


-1) 


3.3849 • 10" 


-6 


3.4058 • 10" 


6 


7.8207 • 10- 


-6 


( 1. 


-1, 


-2) 


2.6184 • 10- 


-7 


2.6498 • 10- 


-7 







( o, 


-3, 


1) 


2.6184 • 10" 


-7 


2.6498 • 10" 


-7 


4.3151 • 10- 


-7 


(-1, 


1, 


-1) 


9.1552 • 10" 


-8 


9.3050 • 10" 


-8 







(-1, 


-1, 


o) 







1.2873 • 10" 


-11 








Table 3: The elements (5 of the core region B, with three distinct estimates ^13, ^14 and 
of the associated densities. The first two were computed using all cycles up to periods 13 
and 14, respectively, and the third using all cycles with denominator not exceeding 200. 
The floating-point numbers represent exact rationals rounded off to 5-digit precision. 
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5 Ordering by denominator 



Let m be a positive integer. We consider the periodic orbits of the map 7 having de- 
nominator m. According to theorem [2] and equation (|19p . these points belong to the set 
H m + B. We are interested in the determination of their period. 

In section [J] we have seen that scaling by Tin maps z^rn into the ring Z[A]/?7iZ[A], which 
for our purpose is a more convenient representation. If a,/3 E Z[A], we denote by (a), 
(a, (3), etc., the ideals they generate. Thus (1) = Z[A]. 

Let £ be a point in Z[A]/(m). The order t(£) of £ is defined as 

= min{£; > 1 : cj fc £ = £ (mod (m))}. (26) 

From theorem [21 we have that every £ G S is the fractional part of at least one periodic 
point; it follows that the order of £ must divide the corresponding period. In fact, from 
conjecture CD it would follow that the period of a periodic point is typically equal to the 
order of its fractional part. Now #H m = m 3 , which gives the crude bound t(£) < m 3 . 

To obtain a sharper bound for t, we must consider the ideal factorization of m in Z[A]. 
We describe it in terms of the factorization of f{x) modulo a prime p. Since Z[A] is the ring 
of all algebraic integers in Q(A), the discriminant of f(x), which is equal to —44, does not 
have any spurious prime divisor, so that the factorization of f{x) modulo p describes the 
ideal factorization of p, without exceptions [24, theorem 27]. The primes that divide the 
discriminant are p = 2, 11; they ramify, that is, f{x) has multiple roots modulo p. For all 
the other primes p, the polynomial f(x) (mod p) has distinct factors. A prime p is inert, 
splits, or splits completely, respectively, if f(x) decomposes into the product of 1,2, or 3 
irreducible factors, respectively. Each possibility occurs infinitely often, and indeed with 
probability 1/3, 1/2, or 1/6, respectively. This follows from Cebotarev's density theorem 
and the fact that the Galois group of f(x) is S3 [28, p. 129]. 

The following result relates the order of a point £ to the ideal factorization of its 
denominator. 

Proposition 4 Let £ G Z[A], and let m be a positive integer. Then t(£) is a divisor ofT, 
where T = T(m) is computed as follows: 

(i) If gcd(n,m) = 1, then T(nm) = hm(T(n),T(m)). 
(ii) If m = p e with p prime and e > 1, we have 
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ramif. 



inert 



T(p e ) 

p e+l {p-l) ifp = 2 
p e (p — 1) if p = 11 

p £ - 1 (p 2 +p + l) ifp = 0,2 (mod 3) 
p^lp 2 +p + l)/3 i/p = l(mod3) 



spirts p 6 " 1 ^ 2 - l)/3 

sp/iis c. 



p e — 1) i/p = 2(mod3) 
p e l (p — l)/3 if p = 1 (mod 3). 



Proof. With reference to the the congruence (|26p . we see that if £ and m are coprime 
(meaning that the ideal identity (£, m) = (1) holds), then upon division by £, we obtain 
lu 1 = 1 (mod (Vn)), independent of £, which shows that is the multiplicative order of 
to in the ring Z[A]/(m). Thus t(£) is a divisor of the order of the multiplicative group of 
this ring. From the decomposition 

ZJA] ZJA] 
(m) ^ J< 

where Jj are the pairwise coprime ideals divisors of (m), we have, in particular, that T(m) 
can be computed as the least common multiple of its value at each summand. 



If £ and m are not coprime, then division by £ leads to a congruence modulo the ideal 
(m)/(£,m). The multiplicative group of the corresponding finite ring is a divisor of that 
considered above, so the same estimates apply. 

Next we consider primary factors m = p e , beginning with the case e = 1. We denote 
by F p fe the field with p k elements, and by F* fe its multiplicative group. We shall use the 
fact that oj = A 3 is a unit of norm 1, and that the prime p = 3 is inert. 



If p is inert, we have the ring isomorphism 

~F 3 
(P) P 

The Galois group of F p 3 over ¥ p is cyclic of order 3, and is generated by the Frobenius 
automorphism T : qh a p '. The orbits of T are the sets of algebraic conjugates. Denoting 
again by to its reduction to F p 3, we find 

^i+p+p 2 = F°(u)ri(u)j#(u) = l. 

The last equality derives from the fact that the terms in the product are the algebraic 
conjugates of uj, and since they are distinct, their product is the norm of oj modulo p. So 
the order of to divides p 2 + p + 1. If p = I (mod 3), then p 2 + p + 1 = (mod 3), and u 
has a cube root in F p 3, which is also a unit of norm 1. In this case the order of uj divides 
(p 2 +p + l)/3. 
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If p splits then 

^i^F p 0F p2 . 

The order of F* divides that of F* 2 , which is equal to p 2 — 1. Because p / 3, we have 
p 2 — 1 = (mod 3), hence u> always has a cube root in F p 2. 

If p splits completely, then 

^y^F p 0F p 0F p . 

The order of the three reductions of u to each finite field is a divisor of the order of F*, 
which is p— 1. If p = 2 (mod 3), there are no further constraints on T(p). If p = 1 (mod 3), 
then <jj has a cube root in F p , and hence the subgroup it generates has index at least 3. 

It remains to deal with the ramified primes 2 and 11. The prime 2 is totally ramified 

f(x) = (x + a) 3 (mod p) p = 2, a = 1 

leading to the ideal factorization (p) = P 3 = (p, A + a) 3 . The ideal P has p 2 incongruent 
points modulo (p), and hence the multiplicative group of Z[A]/(p) has order p 3 — p 2 = 
p 2 (p-l). 

The prime 11 is partially ramified. We find 

f(x) = {x + af(x + b) (mod p) p = 11, a = 3, 6 = 6 

which corresponds to the ideal factorization, (p) = P 2 Pi = (p, A + a) 2 (p, A + b). Hence 

Z[X] 



(P) 



Where R is a ring with p 2 elements, p of which are not invertible. Thus its multiplicative 
group has order p 2 — p = p(p — 1), while that of ¥ p has order p — 1, hence the result. 

Finally, for e > 1, we apply a standard lifting argument. Let k be the largest integer 
for which T(p k ) = T(p). (Such an integer k exists and is effectively computable.) That is, 

oj t (p) = l+p k l3 

where (3 6 Z[A] is such that (f3,p) = (1). Then, using the binomial theorem, we have 

where again (/3',p) = (1). An easy induction on e shows that the order increases regularly 
by a factor p at each step. Thus 

T(p e ) = T{p) p max (0,e-fc)_ 
This completes the proof. D 
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In the totally ramified case, the ^-dynamics is quite regular. Prom the fact that 
(2) = 2i = (2,A + 1) 3 , and u = 1- A- A 2 , we have u - 1 = -A(A + 1) G 2 X . Thus, for any 
£ G Z[A], we have that w£ = £ (mod 2i), and hence every orbit of to modulo (2 e ) consists 
of congruent points modulo 2i. So, if the period is maximal, orbits are cosets of a lattice. 

The theorem states that is a divisor of T, but gives no information as to the actual 
value of T/t. The exact value of is (essentially) determined by the order of the image 
of uj in various finite fields, which cannot be computed in non-polynomial time. One would 
expect that the occurrence of a given ratio T/t will be accompanied by large fluctuations 
and by a limiting probability, which should be maximal for T/t = 1. As a rule, these are 
very hard problems, of the kind mentioned in the introduction in connection with Artin's 
conjecture. 

Proposition 0] generalizes, with appropriate modifications, to any cubic field. What is 
essential is that w be a unit: the fact that our u> is the third power of a fundamental unit 
has merely brought about some specialization in the formulae, expressed via congruences. 
The presesence of the split case in proposition |4] is a consequence of the fact that the 
Galois group of f(x) is the symmetric group. If the Galois group of f(x) is cyclic, the 
split case does not occur, while the other formulae in the proposition remain the same. In 
this case the Cebotarev's densities for the allowed factorizations become 1/3 for primes 
splitting completely, and 2/3 for inert primes. 

We have constructed all periodic points with denominator not exceeding 200. They 
have the form x = £ + (3 with G B and £ G U m <2oo^m> where E m was defined in (fl9j) : we 
obtain 373112717 points in total. Now, every fractional part £ is represented in this set, 
and the total number of fractional parts is given by 



where \i is the Mobius function. (The divisor sum gives the number of lattice points in H m 
that do not belong to any lattice of smallest index.) The difference between the number of 
periodic points and the number of their fractional parts gives the size of the boundary set 
restricted to these denominators (to match the figures, one has to take into account the 
fact that two £-values correspond to three and four (3- values, respectively). In our data 
set the boundary points are one in 10 6 . These findings agree with the corresponding data 
for orbits with increasing period, and underpin our conjecture UJ 
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